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Abstract 



Here we develop a model for the relativistic spring. We examine the effects of revising 
the simple harmonic oscillator to include relativistic momentum and a delayed force 
law. These corrections alter two of the most significant features of the simple harmonic 
oscillator: energy conservation and a constant period independent of initial conditions. 
The relativistic momentum correction, while preserving energy conservation, does not 
have period independent of initial conditions. The delayed force law, while preserving 
period independence, does not conserve energy. Applying both corrections creates a 
solution with increasing amplitude and increasing period, a result that is very different 
from the traditional simple harmonic oscillator. 
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< x > 

Figure 0.1: The simple harmonic oscillator. 

Every freshman physics major is familiar with the simple harmonic oscillator. First 
we introduce Hooke's law in one dimension, F — —kx, where k is the spring constant 
and x is the displacement from equilibrium. Then we invoke Newton's second law, 
mi — p — F, where p is the total time derivative of momentum. From these we 
recover the equation of motion: 

mx = —kx. 



We impose the initial conditions that the mass starts from rest, i.e. x(t)| t =o = 0, and 
starts at initial displacement a, i.e. x(t)| t=0 = a. Using these we solve the ordinary 
differential equation and obtain the position as a function of time, 



(1) 




It is here that we can pause to appreciate two qualities of interest: energy conservation 
and period which is independent of initial conditions. 
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K — ^tnx 2 , so the total energy of the system is 

1,1, 
E = K + V = -kx 2 + -mx 2 . 
2 2 

Plugging in our solution for x(t) we can show that energy is conserved. 

* = + 2 m ( a \/l sin (/I t )) =5 fca2 ' 

which is a constant, so ^ = 0. As there is no change in energy with time, energy is 
conserved. 

Additionally, we find that the period is independent of initial conditions. By inspect- 
ing CO, we note that the simple harmonic oscillator has period given by 

This is only dependent on the parameters of the setup, k and m, and is not dependent 
on the initial displacement a. 

With the simple harmonic oscillator described, we now raise the concerns brought 
up by Einstein's special relativity. Recall that special relativity introduces two postu- 
lates, the principle of invariant light speed and the principle of relativity. While these 
postulates lead to a variety of interesting effects, we will focus on only a few features 
of special relativity which are relevant to our study. We should note that in Einstein's 
theory, energy E is given by 

E — \j p 2 c 2 + m 2 c 4 , 

where p is relativistic momentum, and c is the speed of light. We should also note that 
relativistic momentum is defined as 

mx 

P = j . ■ (3) 

For relativistic momentum we then have that as x — » c, p — > oo, such that exceeding 
the speed of light becomes impossible. So we get a universal "speed limit," c. We also 
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get a new, relativistic form of Newton's second law. From p — F we get 




It is here that we begin to see conflicts with our current notion of the simple harmonic 
oscillator. Our solution CD does not prohibit an oscillator which moves at a velocity 
greater than c. The use of relativistic momentum © will take care of this problem. 

We should note that nothing — not even information — can travel faster than the 
speed of light. The linear restoring force we used earlier, F — —kx, assumes the in- 
stantaneous transfer of information along the spring. We will want to correct this by 
replacing this force with the gradient of a retarded potential, similar to what is done in 
electrodynamics. [1] 

Additionally, special relativity introduces ambiguity in our linear restoring force: it 
is not clear in what reference frame x is determined. Accordingly, we will explore a 
correction to the spring constant k which arises when calculating the force due to the 
displacement in the rest frame of a mass and then transforming into an inertial "lab 
frame." 

Let us proceed to reconcile these issues. 



Chapter 1 

Relativistic Momentum 



We begin our examination of the relativistic spring by introducing relativistic momen- 
tum into Newton's second law. Doing so will ensure that the oscillator does not move 
at a velocity greater than c, our first grievance with the simple harmonic oscillator. 

We will first show that this model conserves energy, then derive a nondimensional- 
ized differential equation and proceed to discuss its numerically determined solutions. 
We will show that this correction gives a period dependent on initial conditions, unlike 
the non-relativistic result. 



We invoke Hooke's law: 



where k is the spring constant, x is the displacement from the origin, and p is relativis- 
tic momentum. Recall that we have relativistic momentum 



1.1 The Model 



mx 



P = 




where m is the mass of the spring, and c is the speed of light. We take then 




(1.1) 
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as our first statement of the model. We will proceed to discuss energy conservation and 
periodicity, and recover from (II. ID a second-order ODE which we will solve. 



1.1.1 Energy Conservation and Periodicity 

In order to arrive at our goal of an energy conservation statement, we wish to write 
this equation as a total time derivative. We will do so by first multiplying both sides by 
x: 




mx 



-kxx . (1.2) 



While the right side falls trivially, 

d (1 _ 

—kxx — > < —kx 

dt [2 

the left hand side requires some coercion. With some perseverance (algebra), we find 
that 



d | mx I d | mc 2 

x- 



dt 



Putting the two results together with (II .2D gives 



d j 1 _ mc 2 

— { -kx 2 + 

dt 2 7772 



Integration by time, where we introduce constant E, yields 



1 „ mc 2 

-kr 2 + — =F. (1.3) 



the desired energy conservation statement. And so, just like the non-relativistic simple 
harmonic oscillator, the model conserves energy. 

Unlike the non-relativistic oscillator however, this model will not have a period 



1.2. Solution 
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independent of initial conditions. If we assume T is independent of a, we can choose 
an initial displacement a such that T < as relativistic momentum requires x < c, it 
is impossible for the oscillator to maintain period T, Thus period must depend on the 
initial conditions. 

We can recover from (11.31) a second order ODE to probe these qualitative observa- 
tions, and we will proceed to do so. 



1.1.2 Getting Our ODE 

Taking the total time derivative of our energy conservation statement (I1.3D gives us a 
second-order ODE, 

mxx 

kxx + : - = 0. (1.4) 

(1 _ (£)2)3/2 

This is not analytically tractable, but can easily be handled using numerical techniques. 
In preparation for numerically solving, we nondimensionalize equation Q1.4D . Let 
t — t T and x — x x, where t and x are dimensionless variables and t and x are 
constants which carry the dimensions time and length respectively. Substituting these 
variables into our second-order ODE gives 



(1.5) 




where it has made sense to let 

/m fm 
- and x = ^-c. 

Here t is merely the fundamental time scale of the problem. Note that t — O) -1 , 
where co is the frequency of oscillation in the non-relativistic SHO problem. From t 
we then get the fundamental length scale of the problem x . From here we can proceed 
to examine the numerical solutions of (I1.5D obtained by Mathematica's Runge-Kutta 
implementation NDSolve. 
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Figure 1.1: Results of numerically integrating (11.51) . We note that the initial conditions 
vary from the classical limit (light gray), to the ultra-relativistic limit (black). Specifi- 
cally: light gray: a = .3x , gray: a = lx , black: a = 1.7x .[2] 



1.2. Solution 
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1.2 Solution 

Let's examine the solutions to (I1.5D with the initial conditions that the mass starts at 
rest with some initial displacement a (i.e. ^ = 0, x(t = 0) = a). We will vary the 
initial displacement a to observe the behavior of the oscillator as we move from the 
classical regime to the ultra-relativistic regime. We see the results of solving numeri- 
cally by Runge-Kutta, courtesy of Wolfram's Mathematica, in figure 11.11 The mass acts 
familiarly in the classical limit-when the initial displacement is significantly less than 
the fundamental length of the problem (i.e. a <C \p^ c )- However as we move into the 
relativistic regime, increasing the initial displacement to greater than the fundamental 
length x , we see a shift in behavior. The velocity function begins to degenerate into 
a square wave as the mass spends more and more time at the maximum speed c. Ac- 
cordingly, the position function begins to look like that of a photon bouncing between 
two mirrors, quickly changing direction when it reaches its maximum displacement, 
and heading back in the opposite direction at speed c. This aligns with our expecta- 
tions, as the substitution of relativistic momentum explicitly imposes a speed limit on 
the model. That is, p — » oo as x — » c. 

We also note that in the classical limit the period depends on m and k, as expected. 
However, as the initial displacement increases, the period increases as well until it is 
entirely dependent on the initial displacement: 

4a 

T — » — . (1.6) 

c 

The period is simply the time it takes at speed c to travel the full wavelength, that is, 
to go from a to —a and back to a. 

We can watch this switch happen in our plot of period versus initial displacement 
in figure [L2l 
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Figure 1.2: Plot of period T by initial displacement a. We see can see the shift from 
the classical solution T = 2n ^j, which is independent of initial conditions, to the 



relativistic solution T = — . 

c 



Chapter 2 
Delayed Forcing 



One of the most central ideas of special relativity is that nothing, not even informa- 
tion, can travel faster than the speed of light. The force law used in the last chapter, 
F — —kx, assumes that the position information travels along the spring instanta- 
neously. We must modify it and delay the force law. We will do so by replacing the 
force with the gradient of a retarded potential, similar to what is done in electromag- 
netism. 

2.1 The Model 



k 




< X\ >< X2 > 

Figure 2.1: Two masses connected by a spring centered about the origin. 

Consider the model from the last chapter. We can think of the fixed point of the 
spring as the source of the spring potential and the mass as a test particle moving in the 
potential. Analogous to the force on one charge due to another in electrodynamics, the 
status of the "source" right now doesn't matter-for the spring potential "news" must 
travel at the speed of light. Instead, the status of the source some time in the past 
informs the potential. However, a fixed source particle (the wall in our prior model) 
subverts this problem by never changing its status, so to see the effects of this we must 
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instead consider a different configuration: two masses connected by a spring. We will 
refer to them as the source mass and the test mass. The news takes a time * s to 

c 

leave the source and reach the test particle, where x t and x s are the locations of the 
test mass and source mass respectively. We can then introduce the retarded time t r 
defined as 

_ l*t(0-*s(£ r )l 

t r t 

c 

By this definition, if news arrives at x t (t), then x s (t r ) is the location of the source mass 
when the news left. Note that we cannot explicitly solve here for t r because this would 
require already knowing the functions x s (t) and x t (t). 

Now we can simply plug t r into the Hooke's law potential to get the retarded po- 
tential, 

y = ifc(x t (t)-x s (t r )) 2 . 

Taking the negative gradient of this potential should give us our force F. In this case, 
we simply take the derivative with respect to x t : 

F = -k(x t (t)-x s (t r )). (2.1) 

Using this force with Newton's second law, we create a set of coupled equations for the 
two masses. First we consider the assignment where, m 1 is the test particle and m 2 is 
the source; then we use m 2 as the test particle and m 1 as the source (figure [27TI) . This 
gives us 

mx 1 = -fc(x 1 (t)-x 2 (t r )), (2.2) 
mx 2 = -k(x 2 (t)-x l (t r )). (2.3) 



With symmetric initial conditions we only require one t r because the two masses are 
at every point indistinguishable, just reflected in space about the origin. Accordingly, 
to simplify our equations we will assume symmetric initial conditions and note that it 
is always possible to find a reference frame that has symmetric initial conditions. We 
choose retarded time 

\x 1 {t)-x 2 {t r )\ 



2.1. The Model 
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Now we have a clear idea of our system, with the ODEs necessary for solving the 
trajectories of the masses. 

At the risk of sounding pedantic, we must make this clear: we are no longer de- 
scribing the behavior of a true-to-life spring, in fact this system is a far cry from your 
everyday slinky; we are simply describing two point masses interacting via a delayed 
spring potential. Accordingly we must check our intuition at the door. We will proceed 
by discussing the qualitative features of this system. 



2.1.1 Energy Conservation and the Retarded Potential 

The system as described by (12. 2D , (12.31) , and (I2.4D appears to not conserve energy by 
construction. Consider a point of maximum displacement between the two masses, 
where we might naively assume we have maximum potential energy and no kinetic 
energy, as in the simple harmonic oscillator. However, in this formulation, at the maxi- 
mum displacement we must have a finite retardation by definition, and so the potential 
is not yet at its maximum. Instead the potential will reach its maximum some finite 
amount of time after the maximum displacement, when the masses have accrued some 
amount of kinetic energy, and so the total energy at that time will be greater than the 
maximum potential energy. Because of this, even though we are deriving the force as 
the gradient of a potential, we cannot think of the force equation as being the time 
derivative of an energy conservation statement. 

2.1.2 Combination of ODEs 

To facilitate discussing the behavior of this system we combine our ODEs into one equa- 
tion for the distance between the masses. In order to do so we introduce displacement 



*(0 = *i(0-* 2 (0- 



(2.5) 



Our equations (CO) . (1231) . and JZ4l) become 



mz{t) = -k{z{t)+z{t r )), 
\\z{t) + z{t r )\ 



(2.6) 



(2.7) 



t, = t 

2 c 
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2.1.3 Nondimensionalization 

Now we need to nondimensionalize these equations. Let t — t T, x 1 — x x 1} and 
x 2 — x x 2 . With these substitutions, our ODEs become 



d 2 x, 
~dz T 

2 



d 2 x 



dT 2 

Here 



-Oi(V) + * 2 Or)) 
-(x 2 (t) + x 1 (t ; .)). 



T r = T-|x 1 (T)-X 2 (T r )|, 

and it has made sense to again assign 



m m 

t ° = \lk and Xo = V¥ c - 

Remember, t is the fundamental time scale of the problem and x is the length scale. 
Similarly, for (I2.7D we let t — t z and z = x z, to get 

d 2 z 

= _ Z ( T ) + Z ( T ) 3 (2.8) 

dz z 

1 

T r = T--|z(T)+ Z (T r )|. (2.9) 



2.1.4 Choosing the Correct T r 

We note that with our definition in (I2.9D . we can have more than one value of T r for a 
given t. The absolute value in our definition yields two possibilities: 

1 1 

T r = t--(z(t) + z(t,.)), and T r = t + -(z(t) + z(t,)). 

From these we get 

11 11 

T r H — z(T r ) = t z(t), and T r z(t,.) = t H — z(t). 



2.2. Solution Method 
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Consider the functions F(t), G(t): 

1 1 

F(t) = t--z(t), and G(t) = t + -z(t). 

For general z(t), F and G may be multivalued; we will show that this only occurs when 
the masses exceed speed c. For single valuedness we have the requirement that F and 
G are monotonically increasing functions (to avoid maxima and minima) . For this we 
need FYt) > and GYt) > 0, which is the case when l^^l < 2. This is true when 
neither mass exceeds the speed of light. Still, we have two possible values for T r : one 
which is less than t and one which is greater. We will choose strictly T r < t as the 
retarded time; the other option is the advanced time and is not physically relevant to 
us. We should note that when the masses exceed the speed of light — and they will 
until we re-impose the correct relativistic form of Newton's second law — we can have 
multiple values with i r < t and no clear idea of which to pick as the retarded time. 
We will see how this can incite some trouble when determining numerical solutions. 



2.2 Solution Method 

To deal with the functional dependence of our ODEs on T r we will implement Verlet 
integration. We will show that basic Verlet integration only requires a slight alteration 
to provide solutions using the retarded potential. [3] 



2.2.1 A Brief Aside on Verlet 

In the Verlet method we begin by discretizing time, taking a grid of size At. If we take 
the Taylor expansion of the position of one of our particles, z(t), at two gridpoints 
t + = t + At and t_ = t — At, we obtain 

z(t)At 2 'z'(t)At 3 

z(t + At) = z(t) + z(t)At + + + 0(At 4 ) 

2 6 

z(t)At 2 'z'(t)At 3 

z(t - At) = z(t) - z(t)At + — + ^(At 4 ). 

2 6 

Adding these together and solving, we can find an expression for z(t + At) 



: z(t + At) = 2z(t) - z(t - At) + z(t)At 2 + <^(At 4 ). 
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And so we get a function for z(t + At) which is only dependent on z(t), z(t — At), 
and z(t). With initial conditions and a description of z we can take steps in time and 
plot the motion of our system. 

2.2.2 Computing T r 

Our function for z(t) requires that we are able to find T r , and now we will describe 
our method for doing so. By definition we have 

1 

T r = T- -|z(T)+z(T r )|, 

which (as mentioned earlier) is not explicitly solvable for T r . But we see that to discern 
the correct x r , we can define the function F 

1 

F(<j) = t - -|z(t) + z(cr)| - a, 

which has the quality that F(o~ — t,.) = 0. In order to find T r we can repeatedly bisect 
the region < a < x r , looking for the point where F(cr) = 0, and taking a at the zero 
crossing of F to be x r . This assumes that there is only one T r < x. In the case that 
the masses are moving faster than the speed of light we are arbitrarily choosing the 
bisection closest to current time x, which may lead to inaccuracies. However, we are 
not really interested in the solution in this case because we know it is unphysical. 

2.3 Solution 

2.3.1 Growth and Periodicity 

Figure 12.21 shows the results of our numerical methods for three different values of 
initial displacement. The solution with the smallest initial displacement is visually in- 
distinguishable from the non-relativistic simple harmonic oscillator, but both of the 
others begin by exhibiting some growth. While we do not have an exhaustive expla- 
nation, we can at least qualitatively motivate the growth by comparing the retarded 
harmonic oscillator to the simple harmonic oscillator. Looking at some arbitrary point 
in the first quarter period of motion some finite time after releasing the mass from rest, 
we can see that the displacement "seen" by the retarded oscillator is greater than that 
seen by the simple oscillator. As this is true for the entire first quarter period, the re- 



2.3. Solution 
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tarded oscillator is accelerated more and is at a greater speed when it makes its first 
zero-crossing. 

At the point of the first zero-crossing, |xiW ~ X2(tr)l _> o and so t r — » t. That is, the 
delay vanishes. Provided x < c, after a zero-crossing the retarded oscillator cannot 
look back beyond the time at which the zero-crossing occurred t zero . To see this we will 
focus our attention to some time after a zero crossing t' — t zero + 5. By definition we 
have the associated retarded time t' r : 

t = t zero + 5 - LlLJ. (2.10) 

c 

We can preserve the ordering t' T > t zem by requiring that 5 > \^LlJ^-il_ Qn the other 
hand we have a lower bound for the time elapsed since last zero-crossing: 5 > rj^- \ . 
Provided > !i^jz^Ol j which is certainly the case if x < c, we have t' > t zero . 

This demonstrates that the retarded time for a given location never goes past the last 
zero crossing unless the mass moves faster than the speed of light. 

At some arbitrary point in the second quarter period, we can see that the retarded 
oscillator "sees" a smaller displacement than the simple oscillator (provided it cannot 
look back past the last zero-crossing, i.e. t' r > t zero ). As this is true for the entire 
second quarter period, the retarded oscillator is decelerated less and thus travels out 
to a greater maximum. For the beginning of the third quarter period, the retarded 
oscillator is looking back beyond the maximum and "seeing" a smaller displacement 
than the simple oscillator, and thus accelerating less. However, it is clear from the 
numerical solutions that this effect is not sufficient to prohibit growth. 

In figure [2721 (c), we can see that once the oscillator begins to break the speed limit 
(x > c), "random lookbacks" begin to occur. That is, the order t' r > t zero is no longer 
preserved, and the retarded oscillator can look back arbitrarily far when determining 
the retarded force. This arbitrary force calculation makes for spontaneous growth and 
decay, which appears as an artifact of the oscillator moving faster than the information 
in the system. The seemingly random growth and decay is the result of demanding a 
speed limit but not enforcing it. 

The important physical result of this section is the loss of energy conservation 
demonstrated by the growth of the oscillations. We can account for this by noting 
that we are not looking at a closed system. By imposing a retarded potential we have 
created the need for a field to mediate the spring force. As we are not accounting for 
the energy and momentum of this field, we have no sense of energy conservation in 
the system. 
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Performing a fast Fourier transform on the numerical solution reveals that periodic- 
ity is independent of initial conditions and constant in time. Comparing the FFTs of the 
data in figure 12.21 displayed in figure 12.31 we see that Fourier transforms have peaks 
at approximately / = V2/2n Hz. This is precisely the frequency that we predict with- 
out any relativistic corrections (where the -/2 comes from coupling the oscillators, see 
equation (12 .6D as T r — » t). This is despite the vastly different initial conditions and 
behavior of the solutions (i.e the growth behaviors shown in figure [2721) . In Chapter 
1 we saw that relativistic momentum changes the periodic behavior without affecting 
energy conservation; now, we see that retarding the potential changes energy conser- 
vation without affecting the period. Naturally, our next endeavor will be combining 
relativistic momentum and the delayed force. 



2.3. Solution 
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Figure 2.2: Results of our Verlet integration for (a) a = 0.001x o , (b) a = 0.050x o , and 
(c) a — 0.01x o ; At = 0.005; 35 000 steps. Trajectories of mass 1 and mass 2 (grey), 
the growth envelope of b (black line). Algorithm in appendix. 
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Figure 2.3: Results of the FFT of our Verlet integration for (a) a = 0.001x o , 
(b) a = 0.050x o , and (c) a = 0.01x o ; At = 0.005; 35 000 steps. Algorithm in ap- 
pendix. 



Chapter 3 
Putting It Together 



In the past two chapters we introduced two corrections independently. First we in- 
troduced the relativistic correction to momentum, and then we formulated a delayed 
Hooke's law force. With the replacement of momentum by relativistic momentum, 
which put a restriction on the maximum speed allowed, we lost the period indepen- 
dence of initial conditions without affecting energy conservation. With our delayed 
force, which restricted the transfer of information in the system, we lost energy con- 
servation without affecting period. Imposing both of these corrections simultaneously 
should produce interesting results. We predict that we will lose both energy conserva- 
tion and simple periodicity. 

We will see the loss of energy conservation for the same reason as in the last chapter. 
Even with the correction to momentum, we are not looking at a closed system, we 
are still neglecting to take into account the field which mediates the force. With the 
addition of the restriction of maximum speed c, we should again see periodicity which 
is not constant. In fact, coupling growth of amplitude with relativistic momentum, we 
should see growth of period over time. 

As we know the growth envelope of amplitude for a retarded potential from Chapter 
2, and we know the period as a function of amplitude for the oscillator in Chapter 1, we 
should be able to predict the behavior of the amalgamation of the models. By plugging 
in the growth envelope of our results from figur e|2.2| into the function in Fig |1.2| we can 
predict the changing period that should result in this final model. In l3.1l we show our 
prediction for period growth as a function of time. 

In addition, it is interesting to include one final relativistic effect. One might take 
the position that the original simple harmonic potential ^kx 2 should be interpreted 
as referring to the length x in the rest frame of the test mass. In this case we get a 
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a 



Figure 3.1: Our prediction from combining figure [2721 and the function in figure [L2l 
See here the period shift that occurs over time due to growth; we should get longer 
and longer period. 



transformation of the spring constant when moving to the inertial "lab frame." 



3.1 The Model 

Combining the delayed force and relativistic momentum, we get the system of equa- 
tions 



dp 1 d I mij 



dt dt l/ttt^ 



= -fc(x 1 (t)-x 2 (t r )), 



and 



dp 2 
dt 



d 
d 



mi, 



-fc(x 2 (t)-x 1 (t r )), 



(3.1) 



3.2. Solution Method 
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where 

1 \x 1 {t) + x 2 {.t r )\ 

t - t __i^ ^J^. (3.2) 

2 c 

3.1.1 Nondimensionalizing 

In order to solve A3. ID and A3. 2D we must nondimensionalize, as usual. Letting t — t T, 
x 1 = x x 1; x 2 = x x 2 , p 1 = p p l3 and p 2 = p p 2 , we have, 

d ' 1 1 = -(*i(0-x 2 (Tr)) 



dp 2 
dz 



= -(x 2 (T)-x 1 (T r )) 3 (3.3) 



where 



T r = t - |x 1 (t)-x 2 (t ) .)|. (3.4) 



We define p 1 and p 2 as 



and p 2 = dT - - (3.5) 



In this case it makes sense to let 

/m fm 
I' x ° = V T c ' and Po = mc ' 

We see the time and length scale of the problem have stayed the same. Additionally, 
we introduce a natural momentum scale p for the problem. 



3.2 Solution Method 

Because of the natural way in which equations (13.31) depend on momenta, it is valuable 
for us to develop a variation of the Verlet method in terms of momenta and positions. 
Here we will describe this. 

For clarity we will walk through the method with reference to one of the masses, 
m 1 . Accordingly, we will make a few notational changes: we will refer to x 1 as x, x 2 
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as y, and p 1 as p. We will represent derivatives with respect to t with the dot notation 
and we will throw out barred notation for nondimensional variables. All together, fore 
example, we will write — * x. 

Our relevant equations are then, from (I3.3D and (I3.5D . 

p = -(x(T)-y(T r )) (3.6) 
x 

(3.7) 



\/l-x 2 

We start with a grid in time with spacing At, as usual, where we define Xj = x(jAt), 
Pj = p(;'At), and t r = rAt. 
We note the definition of p ; , 

Inverting this for x ; and using the midpoint approximation of x ; , we obtain 

X ^~ X ^ = (3.9) 

2At yr^i 

Solving for we get 

2 pAt 

x J+1 = x-_ 1 + (3.10) 

This allows us to find x J+1 given and p.-. To find p ;+1 , we take the midpoint 
approximation of p ; , 

p^ P]+1 ~ Pj -\ (3.11) 
F; 2At 

where by the force equation, we have 

p j = -(x j -y r ). (3.12) 



3.3. Solution 
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Xi(t), x 2 (t) 



t 



Figure 3.2: Qualitative example of combining the effects, with a — 0.5x o , At = 0.001, 
30 000 steps. 

Together, (I37TTT) and (I3TT21) yield: 

p j+1 = Pj _ 1 -2At(x j -y r ). (3.13) 

So with formulae for p ;+1 and x j+1 we can step in time to find our solutions, provided 
we can find r = t r /At. We will find r by again calling a function F, 

F(r) = rAt - jAt + \Xj -x r \, 

and running a bisection routine on F to find F(r < j) = 0. 

3.3 Solution 

The result of this numerical method is shown in figure 13.21 We can immediately see 
the period growth and amplitude growth together, as was expected. As the oscillator 
grows in amplitude, it begins to travel at its maximum speed, and we see the shift from 
sinusoid to triangle wave. So the particles increase their maximum speeds until they 
look like a pair of photons bouncing between two mirrors as the mirrors move apart. In 
figure [331 we compare the period growth to that which we predicted based on a naive 
combination of the two corrections. While qualitatively the results are similar, they do 
have quantitative differences. 





3.4 Rethinking the Spring Constant 



3.4.1 A Modification to our Model 



Consider for a moment how our equations must change as we change our reference 
frame. When we are in the test mass's rest frame we can calculate the total energy E, 

— ? 1 2 

E = mc H — kx . 
2 

Here, ~x is the separation of the two masses in the test mass's reference frame and k is 
the spring constant in this frame. When we change to the lab frame (or any inertial 
reference frame), we have E — , To recover our formula for total energy in the 
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lab frame, we must have 



This gives 



where 




That is, the spring constant k in the lab frame is different from the spring constant k in 
the test mass's frame. If we assume k is constant, this implies another modification to 
the force equation. With a corrected spring constant plugged into our solution, we can 
see that the growth in amplitude is slowed down (figures 3.4 and 3.5). However, the 
behavior is nearly identical. 
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a 



Figure 3.4: Period shift from solution before (black) and after spring constant correc- 
tion (black dashed). 




Figure 3.5: Qualitative example of the correction applied to figure [3T2l before correc- 
tion (gray), after correction (black). Where, a = 0.5x o , At — 0.001, 30 000 steps. 



Closing Remarks 



The simple harmonic oscillator undergoes a variety of relativistic corrections which 
alter its most familiar features such as energy conservation and a period independent 
of initial conditions. 

In Chapter 1 we implemented relativistic momentum to insure that the oscillator 
never moved faster than the speed of light. While this correction preserved energy 
conservation, it created a period dependent on amplitude. Imposing the "speed limit" 
while increasing the amplitude made the mass's trajectory look increasingly like that of 
a photon bouncing between two fixed mirrors. 

In Chapter 2 we introduced a delayed force in order to avoid having information 
move faster than the speed of light. While this correction had no affect on the period, it 
did cost us energy conservation. This delay implied that information propagate along 
the spring, but we weren't accounting for the energy of the field which mediated the 
force, and so we saw growth. 

In Chapter 3 we found that the amalgamation of these two effects gave us a predic- 
tion of how the combined model would act — an oscillator with increasing amplitude 
and period. Our prediction qualitatively lined up with what we observed numerically, 
although there were small differences. Finally we considered allowing the spring con- 
stant to transform with reference frame shift. This produced only slightly different 
numerical results. 

To continue this project it would be of value to continue stress testing the numerics 
of Chapter 3 and searching for the reason the prediction and numerics of chapter 3 do 
not perfectly align. As an extension, it may be advisable to formulate the field which 
governs the force described in Chapter 2. By doing so, and taking the dynamics of the 
field into account, it may be possible to recover energy conservation. 



Appendix A 
Mathematica Code 



A.1 Bisect While 

This is the Code that is used to find the value of T r : 

BisectWhile [F_, maxz_, eps_] : = 

Module [{xm, xl , xr, vam, val, var, retval, i, pos, k}, 
k = 0; 

(*If [maxz==l, Return [1]] ;*) 
retval = -1; 
xl = 1; 
xr = maxz ; 

xm = Floor [(xl + xr)/2]; 
val = F [xl] ; 
var = F [xr] ; 
vam = F [xm] ; 

If [Abs [F [xl] ] > Abs[F[maxz]] , pos = maxz + 1, 
pos = 1] ; Ochoose the smallest of the past or the current time*) 

If [val var > 0, retval = pos,(*if no zero crossings, 
choose the smallest of the beginning or the current time*) 
While [Abs [vam] > eps && k < 10000, 

If [val vam < 0, xr = xm; , xl = xm;] ; 

xm = Floor [(xl + xr)/2]; 

val = F [xl] ; 

var = F [xr] ; 

vam = F [xm] ; 
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retval = xm; 

]; 

k = k + 1; 
]; 

Ofirst two fix the problem introduced by the floor function, 
i.e. look left and right . the third moves the retval to the 

head if it has been flagged with 
If [retval == -1 I retval == xr, retval = retval, 

If [Abs[F [retval + 1]] < Abs [F [retval] ] , retval = retval + 1] ] ; 
If [retval == -1 I retval == xr, retval = retval, 

If [Abs [F [retval - 1]] < Abs [F [retval] ] , retval = retval - 1]]; 
If [retval == -1, retval = pos] ; 
Return [retval] ; 

]; 



A.2 Verlet 

This is the Code that is used to implement Verlet with the delayed force : 

Verlet [xlnitial_, deltaT_, N_] := 
Module [{tol, i, tr, posl, pos2, z, max, xl, x2, fVals}, 
max = {{0, xlnitial}}; 
z = 0; 

tol = 10 deltaT; 
xl = {xlnitial, xlnitial}; 
x2 = {-xlnitial, -xlnitial}; 
For[i = 2, i<=N, i = i + l, 
fVals[j_] : = 

j deltaT - Length [xl] deltaT + Abs [xl [ [Length [xl] ] ] - x2[[j]]]; 
tr = BisectWhile[fVals, Length [xl] - 1, 10 deltaT]; 
posl = (2 xl[[i]] - 

xl[[i - 1]] - (deltaT~2) (xl[[i]] - x2[[tr]])); 
pos2 = (2 x2[[i]] - 

x2[[i - 1]] - (deltaT~2) (x2[[i]] - xl[[tr]])); 
xl = Join[xl, {posl}]; 
x2 = Join[x2, {pos2}] ; 



A.3. Momentum Verlet 
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If[(xl[[i - 1]] - xl[[i - 2]])*(xl[[i]] - xl[[i - 1]]) < 0, 
max = Join [max, {{deltaT*(i - 1), Abs[xl[[i - 1]]]}}]]]; 
Return [{{xl , x2} , max}]] 

A. 3 Momentum Verlet 

This is the Code that is used to implement Verlet with the delayed force and Relativistic 
Momentum: 

ComboVerlet [xlnitial_, deltaT_, N_] := 

Module[{tol, i, posl, pos2, z, zero, xl, x2, pi, p2, k, c, m, moml, 
mom2, j, fVals, tr} , 
xl = {xlnitial, xlnitial}; 
x2 = {-xlnitial, -xlnitial}; 
pi = {0, 0}; 
p2 = {0, 0}; 
zero = {}; 

For[i =2, i<=N, i=i+l, 
fVals[j_] : = 

j deltaT - Length [xl] deltaT + Abs [xl [ [Length [xl] ] ] - x2[[j]]]; 
tr = BisectWhile [fVals, Length [xl] - 1, 10 deltaT]; 
moml = pi [[i - 1]] - 2 deltaT (xl[[i]] - x2[[tr]]); 
mom2 = p2[[i - 1]] - 2 deltaT (x2[[i]] - xl[[tr]]); 
posl = (2 pl[[i]] deltaT) /Sqrt[l + pl[[i]]~2] + xl[[i - 1]]; 
pos2 = (2 p2[[i]] deltaT) /Sqrt[l + p2[[i]]~2] + x2[[i - 1]]; 
xl = Join[xl, {posl}]; 
x2 = Join[x2, {pos2}] ; 
pi = Join [pi, {moml}]; 
p2 = Join[p2, {mom2}] ; 

If [xl [ [i] ] *xl [ [i - 2]] < 0, zero = Join [zero, {deltaT* (i - 1)}]]; 

]; 

Return [{{xl, x2} , zero}]; 

]; 
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